Cavitation and bubble collapse in hot asymmetric nuclear matter. 
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Abstract 

The dynamics of embryonic bubbles in overheated, viscous and non-Markovian nuclear matter 
is studied. We show that the memory and the Fermi surface distortions significantly affect the 
hinderance of bubble collapse and determine a characteristic oscillations of the bubble radius. 
These oscillations occur due to the additional elastic force induced by the memory integral. 
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1. Growing or collapsing bubbles in an overheated (undercompressed) fluid has been a 
intriguing problems of physics over a long period of time jjj. In nuclear physics, the bubbles 
may be formed in the expanded nuclear matter in a heavy-ion reaction then providing the 
break-up of the matter into fragments 

esq. 

The formation of embryonic bubbles in 
hot nuclear matter due to the quantum and statistical fluctuations has been studied earlier 
in Refs. The goal of the present paper is the analysis of the dynamic 

evolution and the collapse of the bubble in an overheated nuclear matter. 

In a classical liquid, the dynamics of the bubble of radius R is derived mainly from 
the thermodynamic potential <&(R), which reaches a maximum value for the critical radius 
R = R* [l2|. More complicated bubble dynamics occurs in the case of a Fermi liquid 
where the dynamic distortion of the Fermi surface produces an additional pressure tensor 
13 , Q|. Moreover, the Fermi surface distortion and the interparticle collisions lead here 



to the non-Markovian equations of motion for the relevant collective variables jig . Below 
the non-Markovian dynamics will be applied to the problem of the collapse of bubbles. The 
non-Markovian equation of motion for the bubble radius R will be derived in a hot nuclear 
Fermi-liquid. However, the final macroscopic equations of motion including the memory 
kernel can be applied to a wide number of non-Newtonian liquids. 

2. We will consider the dynamics of the bubble (cavitation) with an arbitrary undercritical 
size R < R*, where the critical radius R* of the bubble in an overheated liquid at the 
temperature T. The critical radius R* depends on the overheating temperature AT = T—Tq, 
where To is the boiling temperature in the case of plane geometry. For fixed values of 
temperature T and liquid phase pressure Pn q the critical radius R* is derived from the 
condition of thermodynamic equilibrium of the vapor bubble with the surrounding liquid 



and it is given by 
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R* = _^5_. (i) 
p vap 0AT 



Here, a is the surface tension coefficient, p vap is the particle density in the bubble and <fi is 
the latent heat of evaporation. Considering the dynamics of the surrounding Fermi liquid, 
one can reduce the collisional kinetic equation for the phase-space distribution function to 
Euler-like equations for the velocity field u(r, t) in the form (for details, see Refs. Q,Q,Q) 



m^-u v + mpiu^^) u v + -V V V + V„-^ = — V M PL + T„, ext , (2) 
at p dp p ' 
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where p = p(r,t) is the particle density, e po t = e po t(r, i) is the potential energy density and 

V = V(r, t) is the pressure caused by both the thermal and Fermi motions of nucleons. The 
pressure tensor P' = V' v Jr,t) in Eq. (J2J) is caused by the Fermi surface distortion. In the 
case of the most important quadrupole distortion of the Fermi surface, the pressure tensor 

V satisfies the following equation jis| 

J^ + ^J^ = -^ with A ^ = V -Xm + V M x, - |^VaXa, (3) 
where r is the relaxation time due to interparticle collisions and \ v = x„(r,t) is the dis- 
placement field which is related to the velocity field by uJr,t) = dxu{ r ,t)/dt. In the case of 
a spherical bubble, the displacement field is given by [?], [3| 



Xu{r,t) = R 3 r„/3r 3 for r > R. (4) 

We point out the quadrupole distortion of the Fermi surface does not distort the spherical 
shape of the bubble if the displacement field Xu( r ,t) is given by Eq. (J5J). 

The external force F ext in Eq. (J2j) is caused by the vapor pressure on the liquid. The 
external force on the bubble surface in the radially outward direction per unit area is given 



by naii 

Fext = -Pvap ~ Pqi (5) 

where P vap is the vapor pressure within the bubble and P is the pressure of the saturated 
vapor with respect to a plane surface. Note that the condition of equilibrium between the 
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bubble vapor and the liquid is expressed by P — Pvap and F ext = 0, see Ref. 

The Euler-like equations (j2J) can be reduced to the equation of motion for the radius R(t) 
of the bubble. Multiplying Eq. (J2J) by x^(r,t) = R 2 r u /r 3 , summing over u, integrating over 



r-space and using Eqs. (JSJ), we obtain the non-Markovian equation for the variable R(t) 
BR + 

where 



BR+ 2dR R + I{R] t) = ~dR~ (1 " P ^/ P °)' 



I(R;t) = B [ dt'R{t') exp[(t' -t)/r] K,{t,t') (7) 
is the memory integral. The inertia! parameter B = B(R) and the memory kernel K,(t,t') 



in Eqs. © and (J7J) are given by p, [Tfl, |20 ] 
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FIG. 1: Dependence of the thermodynamical potential A<I> (see Eq. (|10[) ) of metastable (over- 
heated) asymmetric nuclear matter on the radius of the bubble. 

where po is the bulk density of the nuclear matter and ep is the Fermi energy. Denoting the 
value of thermodynamic potential $ corresponding to the absence of the bubble by $o> we 
have 

$ = $(i?) = $ + A$(R), (9) 

where A$(i?) is the change of the thermodynamic potential due to the variation of the 
bubble radius. The value of A$(i?) is given by |3] 

A*(J2) = 4vra (tf - ^ . (10) 

The thermodynamic potential $(-R) for R = R* does not have a minimum value, as in the 
case of equilibrium, but a maximum value, see Fig. 1. The thermodynamic potential of Eq. 
(JTUJ) corresponds to the metastable liquid phase with /ii iq > /i vap , where /zi iq and /i vap are the 
chemical potentials for the liquid and the bubble vapor, respectively. 

Note that in the Markovian limit at r — > the equation of motion (JH1) is transformed to 
the generalized Rayleigh-Plesset equation of classical bubble dynamics 19]. 
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The critical radius R* in Eq. (fTUjl depends on the temperature T of liquid-gas phase 
transition, the overheating temperature AT and the isotopic asymmetry parameter Xn q 
(X = (p n — Pp) I (p n + Pp) , where p n and p p are the neutron and proton densities, respectively), 
sec Eq. © and Refs. Q,Q- 

To evaluate the critical radius R* we point out that the 



caloric curve has the plateau region at the temperature T 



2l\ . Below we will adopt T = 4 



MeV and Xn q = 0.1. If one assumes the process of isobaric heating for the description of the 
plateau region in the caloric curve at To = 4 MeV, the order of magnitude of the pressure 
Pq should be Pq ~ 10~ 3 MeV/fm 3 for this process [2jj. Under these conditions, one obtains 
from Eq. JQ iJ* = <?/AT with q ph 4.1 MeV • fm for a symmetric nuclear matter. The 
parameter q grows slightly with the asymmetry of the nuclear matter. 

In the case of boiling of asymmetric nuclear matter, the vapor asymmetry X v? p signifi- 



cantly exceeds the corresponding liquid asymmetry Xi iq (at Xi iq > 0), see Ref. 21]. This is 
a feature of the nuclear matter where the structure of the isospin symmetry energy provides 
the condition \p n \ < \p p \ and the preferable evaporation of neutrons. Using this fact, we 
will neglect the proton fraction in the bubble vapor. The presence of the noncondensable 
neutron vapor within the bubble then gives a contribution to the term ~ (1 — P vap /Po) in 
Eq. ©. (Note that the analogous term in the classical Rayleigh-Plesset equation is due 
to the assumption that the bubble contains some quantity of contaminant noncondensable 
gas 19].) The vapor pressure in the bubble includes both quantum and thermal contribu- 
tions. Since the temperature region of interest is T e^, the thermal pressure is relatively 
small and for the degenerate neutron gas one can use the Thomas-Fermi approximation with 
Tvap/To = (R*/Rf inEq. ©. 

The surface tension coefficient a in Eq. (fTUj) is temperature dependent. We will use the 
following expression for a = cr(T) 
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a(T) = ff (0) [(T c 2 rit - T 2 )/(T c 2 rit + T 2 )] 5/4 , (11) 
where T crit = 14.6 MeV is the critical temperature for infinite nuclear Fermi-liquid, associated 



with the SKM interaction 



2l|. 



3. Below we will consider the collapse phase, i.e. the descent from the top of the barrier in 
Fig. 1 toward R < R*. In general, an accurate evaluation of the bubble dynamics requires to 
include a thermal boundary condition at the bubble wall to provide the possible changes of 
the pressure P vap of the bubble vapor and its temperature (l9|. We will neglect the thermal 



effects, assuming that the main contribution to the pressure of the bubble vapor is due to 
the quantum motion of nucleons because of T <C ep. If the collapse process is fast enough 
and the condensation of the neutron's vapor within the bubble is negligible, the force F ext 
of Eq. (j5J) has then to be taken into consideration in Eq. (j2J). In an opposite case of a slow 
collapse process, the possible growth of the pressure P vap at the decrease of the bubble radius 
R(t) is compensated by the condensation of the vapor providing that P vap = P = const and 
F ext = 0. 

The non-Markovian equation © can be solved if one rewrites it as a set of two connected 
equations for the bubble radius R(t) and the memory integral I(R; t) with the relevant initial 
conditions. The bubble starts to collapse from the metastable state (point A on the top of 
barrier in Fig. 1) at -R(O) = R*. Taking the overheating temperature AT = 0.5 MeV, we 
obtain from Eq. (0) i?(0) = R* = 8.2 fm. This critical radius R* = 8.2 fm seems too large 
for the finite nuclei. We point out, however, that the inclusion of the Coulomb forces, which 
is ignored in our consideration of the infinite nuclear matter, decreases the value of R* (f| Q| • 
Note, however, that our consideration cannot be extended to very small values of R* because 
for a small enough bubble the finite size effects, in particular the finite diffuse layer for the 
vapor density inside the bubble, can be important 0. Using the above mentioned initial 
condition for R(0), we restrict ourselves to an infinite nuclear matter and expect that the 
main conclusions concerning the Fermi motion effects on the collapse of the bubble in hot 
Fermi liquid can be also applied to finite nuclei. 

The initial velocity R(0) can be derived using the initial kinetic energy -Ekin,o- Assuming 
the equipartition of energy over degrees of freedom at R = R*, we use -Ekin.o = (^o + AT)/ 2 
and obtain R(0) = — y/2E\a n t a /B{R*) = —2.9 • 10 _3 c, where sign "— " is used because of the 
collapse process with R(t) < R*. 

In Fig. 2, we show the time dependence of the radius R(t) in the presence of the noncon- 
densable vapor of neutrons with F ext ^ (solid lines). Here we have adopted ep = 37 MeV, 
cr(0) = 1.1 MeV • fm -2 and po = 0.16 fm -3 . For the relaxation time r we have used r = tq/T 2 
with To = 850 fm/c |23J]. As seen from Fig. 2, the collapse of the bubble is accompanied by 
oscillatory behavior of its radius R(t). These oscillations occur due to the memory integral 
in Eq. (jUJ) and disappear in the limit of Markovian motion (see dashed line in Fig. 2). 
As seen from Fig. 2, the collapse process stops at a certain values of the bubble radius, 
-Rmin = 0.28 R*. The presence of non-zero i? m i n ^ is due to the pressure of the noncon- 
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FIG. 2: Dependence of the radius R(t) of the collapsing bubble on the time obtained for two 
regimes: non-Markovian motion, Eq. 0, (solid line) and Markovian (no memory) motion, Eq. 
()13|) . with the friction coefficient 7 from Eq. Q14|) (dashed line). The inset shows the collapse 
process near the top of barrier of &(R). The values of To, AT and X\[ q as in Fig. 1. 

densable neutron vapor in the bubble. The value of -R m i n is derived by the condition of the 
compensation of the adiabatic forces in Eq. (JBJ at t ^> r. 

The memory integral in Eq. © leads to a significant delay in the collapse process. To 
show this effect we will derive the friction coefficient 7 from the memory integral I(R; t). In 
a immediate vicinity of the critical radius R*, the solution to Eq. (JHJ) takes the form of a 
superposition of exponential functions exp(Ajt) (i = 1, 2 and 3) with eigenvalues \ obtained 
as solutions to the secular equation 



Here k = -4vr(2cr - 5P R*), k = (32/5)7rp e F R* and the mass coefficient B* = B(R*) is 
taken from Eq. (|8j) at R — R*. In the case of the zero-relaxation-time limit, r — > 0, one 



by the static stiffness coefficients k. In the opposite case of rare collisions, r —>■ 00, the 

7 




(12) 



obtains from Eq. ()12|) the motion with A = ±y[ 




solution to Eq. (|T2"j) leads to a motion with A = ±za/(— \k\ + k)/B*, where the additional 
contribution to the stiffness coefficient, k, appears because of Fermi surface distortion effect. 
In both limits r — ► and r — > oo, the macroscopic equation of motion (JHJ) is reduced to the 
usual Markovian (no memory) equation with a friction coefficient 7 

mp RR + ^mp R + fy R=-^(i-^j- P h _ j . (13) 

The friction coefficient 7 is related to the relaxation time r as [15] 

^^ b "ttww- (14) 

where ujf = is the characteristic frequency for the eigenvibrations caused by the 

Fermi surface distortion. The corresponding solution for R(t) obtained from the Markovian 
equation of motion (|13j) with friction coefficient 7, from Eq. (|14jl . is shown in Fig. 2 by the 
dashed lines. As seen from Fig. 2 the presence of memory effect (solid line) strongly hinders 
the collapse process. If the collapse process becomes slow enough the growth of the vapor 
pressure P vap is compensated by the condensation of the vapor and the contribution from 
F ext is negligible in Eq. (JHJ). The collapse process then leads to the disappearance of the 
bubble. 

In Fig. 3 we have plotted the dependence of the time, t c i ps , required for the total collapse 
from R = R* to R = 0, on the relaxation time r, for the bubble in the absence of the force 
F ext (i.e., neglecting the term ~ P in Eqs. (JBJ) and (fT3j) ). 

In Fig. 3 we compare the results for t c \ ps for both the non-Markovian motion derived by 
Eq. © (solid line) and the usual Markovian (no memory) one given by Eq. (|T3*j) (dashed 
line). As seen from Fig. 3, the behavior of t c i ps is changed significantly due to the memory 
effects. A delay in the collapse of the bubble occurs due the non-Markovian (memory) effect 
for the large values of the relaxation time. 

In the limit of a non- viscous liquid, r — >• 0, we obtain for both motions t c i ps ,o = 4 • 10 -22 
s. This result can be compared with the classical one of Rayleigh t c ips,R [24]. To do that, we 
note that under the condition used above for T = 4 MeV, P ~ 10~ 3 MeV/fm 3 , X liq = 0.1 
and R* = 8.2 fm, the Fermi liquid is undercompressed by a pressure AP w 0.25 MeV/fm 3 . 
The Rayleigh collapsing time t c i ps ,R is given by |3] 



idp „ R = 0.915 (^>) 



1/2 
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FIG. 3: The collapse time t c i ps for the metastable bubble in an overheated Fermi liquid as a 
function of relaxation time r for different temperatures To = 4 and 6 MeV shown near the curves 
at fixed asymmetry parameter Xn q = 0.1. The solid and dashed lines represent the calculations 
with non-Markovian and Markovian equations of motion respectively. The values of AT and X\[ q 
are the same as in Fig. 1. 

with the yield t c i pSi R = 6.1 ■ 10~ 22 s which is in a good agreement with t c ips,o- 

4. In conclusion we note that the collapse of the bubble in an overheated (undercom- 
pressed) Fermi liquid is strongly influenced by memory effects, if the relaxation time r is 
large enough at ojp r > 1. In this case, the collapse of the bubble is accompanied by char- 
acteristic shape oscillations of the bubble radius R(t) (see Fig. 2) which depend on the 
memory kernel JC(t,t') and the relaxation time r. These oscillations appear due to the non- 
adiabatic elastic force induced by the memory integral. The non-adiabatic elastic force acts 
against the adiabatic force —d§/dR (see Eq. (|0J)) and hinders the collapse of the bubble. 
In contrast to the case of the Markovian motion, the delay in the collapse is caused here by 
both conservative elastic and friction forces. 

We should like to point out that the study of the cavitation in a two-component nuclear 
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matter could represent a stimulating area of interest due to the fact that the vapor inside 
embryonic bubbles is strongly asymmetric system with (N — Z)/ (N + Z) w 1, where N and 
Z are the number of neutrons and protons, respectively. The collapse of such extremely 
neutron rich system in which the energy is focused in a small region could show up new 
phenomena in the clusterization of nucleons. 

The main goal of this paper was to show the existence of a new phenomena of shape 
oscillations of the collapsing bubbles in a hot nuclear Fermi liquid. This phenomena does 
not exist in classical (Newtonian) liquids. The shape oscillations of the collapsing bubble 
could induce a quite stimulating search of the accompanied characteristic 7-quanta emission. 
That provides, in principle, the possibility for the measurement of the temperature of the 
first kind phase transition through the measurement of the energy and the damping of the 
corresponding resonances in the 7-quanta spectrum. 

Our consideration was restricted to a nuclear matter for maximum simplicity. Both 
effects of Coulomb forces and nuclear surface tension must be taken into account in the 
case of bubble dynamics in finite nuclei. The presence of Coulomb forces enhances the 
bubble formation and decreases the critical radius R*. Note, however, that the minimum 



value of the radius R*, for which our approach is applicable, is determined 
a/R* <C 1, where a = 0.5 -j- 1 fm is the surface thickness of the bubble 
dynamics in hot finite nucleus requires further investigations. 



3y the condition 
. The bubble 



The author thanks H. Hofmann and R. Hilton for valuable suggestions and the Physics De- 
partment of the Technische Universitat Miinchen for the nice hospitality. Financial support 
by the Deutsche Forschungsgemeinschaft under contract 436 UKR 113/66/0-1 is gratefully 
acknowledged. 



[1] D. Lohse, Physics Today 56 (2003) 36. 

[2] J.P Bondorf et al., Nucl. Phys. A 443 (1985) 321; A 444 (1985) 460. 
[3] H.W. Barz et al., Nucl. Phys. A 448 (1986) 753. 
[4] K. Sneppen, Nucl. Phys. A 470 (1987) 213. 

[5] H. Schulz, D.N. Voskresensky and J.P. Bondorf, Phys. Lett. B 113 (1983) 141. 

10 



[6] A.H. Blin et al., Phys. Lett. B 182 (1986) 239. 

[7] A.H. Blin, B. Hiller, H. Reinhardt and P. Schuck, Nucl. Phys. A 484 (1988) 295. 

[8] J.P. Bondorf, R. Donangelo and K. Sneppen, Phys. Lett. B 214 (1988) 321. 

[9] R. Donangelo, K. Sneppen and J.P. Bondorf, Phys. Lett. B 219 (1989) 165. 

[10] A.J. Santiago and K.C. Chung, Nouvo Cimento A 105 (1992) 1729. 

[11] L.D. Landau and E.M. Lifshitz, Statistical Physics (Pergamon, Oxford, 1958). 

[12] J. Frenkel, Kinetic theory of liquids (Clarendon Press, Oxford, 1946). 

[13] A. A. Abrikosov and I.M. Khalatnikov, Rep. Progr. Phys. 22 (1959) 329. 

[14] V.M. Kolomietz and H.H.K. Tang, Phys. Scripta 24 (1981) 915. 

[15] V.M. Kolomietz, S.V. Radionov and S. Shlomo, Phys. Rev. C 64 (2001) 054302. 

[16] V.M. Kolomietz and S. Shlomo, Phys. Rep. 690 (2004) 133. 

[17] H. Lamb, Hydrodynamics, Art. 91a (Dover Publications, New York, 1945). 

[18] M.S. Plesset, ASME J. Appl. Mech. 16 (1949) 228. 

[19] C.E. Brennen, Cavitation and bubble dynamics (Oxford University Press, Oxford, 1995). 

[20] V.M. Kolomietz, A.I. Sanzhur and S. Shlomo, Phys. Rev. C 68 (2003) 014614. 

[21] V.M. Kolomietz, A.I. Sanzhur, S. Shlomo and S.A. Firin, Phys. Rev. C 64 (2001) 024315. 

[22] D.G. Ravenhall, C.J. Pethick and J.M. Lattimer, Nucl. Phys. A 407 (1983) 571; 

D.G. Ravenhall, C.J. Pethick and J.R. Wilson, Phys. Rev. Lett. 50 (1983) 2066. 

[23] P. Danielewicz, Phys. Lett. B 146 (1984) 168. 

[24] Lord Rayleigh, Phil. Mag. 34 (1917) 94. 



11 



